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Abstract 

We study the evolution and scaling of the entanglement entropy after two types of quenches for 
a 2+1 field theory, using a conjectured holographic technique for its computation. We study a 
thermal quench, dual to the addition of a shell of uncharged matter to four dimensional Anti-de 
Sitter (AdS4) spacetime, and study the subsequent formation of a Schwarzschild black hole. We 
also study an electromagnetic quench, dual to the addition of a shell of charged sources to AdS4, 
following the subsequent formation of an extremal dyonic black hole. In these backgrounds we 
consider the entanglement entropy of two types of geometries, the infinite strip and the round 
disc, and find distinct behavior for each. Some of our findings naturally supply results analogous 
to observations made in the literature for lower dimensions, but we also uncover several new 
phenomena, such as (in some cases) a discontinuity in the time derivative of the entanglement 
entropy as it nears saturation, and for the electromagnetic quench, a logarithmic growth in 
the entanglement entropy with time for both the disc and strip, before settling to saturation. 
We briefly discuss the possible origin of the new phenomena in terms of the features of the 
conjectured dual field theory 
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1 Introduction 



The entanglement entropy of a quantum system has emerged as a useful probe in a number of 
fields. For example, the condensed matter community have recognized its utility as a diagnostic of 
phase transitions, especially in contexts where an analogue of a Landau-Ginsburg order parameter 
is either not accessible or available. The quantum information community employ it as a tool for 
assessing such things as the computational complexity of a system (see ref. [I] for a review). The 
quantum gravity community have used it to address important issues such as the entropy associated 
with a region of spacetime such as a black hole, where (for example) the important relation between 
the entropy of an excised region of spacetime (such as the interior of a black hole) and the area of 
the spacetime surrounding it can be motivated in terms of entanglement entropy (see e.g. ref. [2].) 
In such studies, it is very important to be able to extract information about how the entanglement 
entropy scales with the size of a given system, or evolves as a function of time. 

Given a system, consider a region or subsystem which we can call A, with the remaining 
part of the system denoted by B. A definition of the entanglement entropy of A with B is given by: 

Sa = -TtU (PA^PA) , (1) 

where pj± is the reduced density matrix of A given by tracing over the degrees of freedom of £>, 
p_4 = Trg(p), where p is the density matrix of the system. When the system is in a pure state, 
i.e., p = the entanglement entropy is a measure of the entanglement between the degrees 

of freedom in A with those in B. If the system starts in a pure state |^(0)) and the evolution 
of the system is unitary, then the system at subsequent times is also a pure state \^>(t)), and so 
the evolution of the entanglement entropy describes the evolution of quantum entanglement in the 
system. 

The evolution of entanglement entropy in a (l+l)-dimensional field theory after a quench 
was studied in ref. [3]. The system for t < is in the ground state (with a mass gap), and at 
t = 0, the magnetic field is set to the critical point of the theory. Dividing the system into A (a 
line segment of length i) and B (the rest), the analytical results for the system showed that the 
entanglement entropy grows linearly with time until it saturates at approximately t = Ijl. Ref. [3] 
explained this as follows: at t = 0, the quench results in excitations that propagate to the left and 
to the right. Since the system before the quench had a mass gap, only excitations produced at 
around the same local region will be entangled. Contributions to the entanglement entropy occur 
when a left /right mover remains in region A whereas the right /left mover is in region B. So since 
at time t, the left-right pair are separated by a distance 2t, saturation is reached when excitations 
produced at the center of A reach B which should be at 2t ~ t. It is of interest to study more such 
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systems and also to go beyond two dimensions. It is hard to make progress for many interacting 
field theories of interest, and so a search for powerful techniques from as broad a range of sources 
seems prudent. One such source of tools is the AdS/CFT correspondence [U EJ E]. 

A prescription for calculating the entanglement entropy [7] and its evolution [8] in the 
context of AdS / CFT provides a new way to calculate the entanglement entropy using geometrical 
techniques that are classical in spirit (for a review see ref. [9]). In an asymptotically Anti-de Sitter 
(AdS) geometry, consider a slice at constant AdS radial coordinate z = a. Recall that this defines 
the dual field theory (with one dimension fewer) as essentially residing on that slice in, in the 
presence of a UV cutoff set by the position of the slice. Sending the slice to the AdS boundary at 
infinity removes the cutoff (see ref. |10j for a review). On our z = a slice, consider a region A. Now 
find the minimal-area surface 7.4 bounded by the perimeter of A and that extends into the bulk 
of the geometry. Then the entanglement entropy of region A with B is given by: 

Area(7^) 

where Gn is Newton's constant in the dual gravity theory. This prescription for the entropy 
coincides nicely with (l+l)-dimensional computations of the entanglement entropy, and has a 
natural generalization to higher dimensional theories. Note that there is no formal derivation of the 
prescription. Steps have been made, such as in refs. I12j. but they are not complete. However, 
there is lots of evidence for the proposal. See e.g., refs. [T3J EH EI3 ESI E3 EE]- A review of several 
of the issues can be found in ref. [9] . Further progress has been made recently in ref. |19| . Although 
there is no derivation available, in this paper we shall assume that this holographic prescription 
does give the correct result for the entanglement entropy in higher dimensions. Our results are 
interpreted in line with this. In particular, we find features that conform to expectations, and in 
addition we find new phenomena that we regard as predictions for the behavior of the field theory 
entanglement entropy. 

Interesting recent papers, ref.jEl E0], presented a study of the time evolution of a (1+1)- 
dimensional system after a thermal quench, by working in the AdS3 geometry dual to it. Since, 
via the AdS/CFT correspondence, a thermal state in the field theory is dual to a black hole on 
the gravity side[2TJ, the authors started with AdS3, perturbed it with a shell of matter, and then 
followed the subsequent formation of a black hole that took place, computing the entanglement 
entropy as a function of time using the above geometric prescription. They used the exact AdS3 
Vaidya metric to describe the formation process. 

Their results confirmed a number of key features seen in the CFT computation of ref. [3]. In 
particular, the linear rise of the entanglement entropy with time, and the scaling of the saturation 
time with the size of the system showing that the effective propagation speed of the entanglement 
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was the speed of light. (Among the many interesting features of such computations is the fact that 
the minimal surfaces involved must probe the region behind the apparent black hole horizon.) 

The goal of the work we report on was primarily two-fold. First, we wished to examine 
the question of higher dimensions, and so we used an AdS4 setting, relevant to the study of dual 
2+1 dimensional field theory. The question of how the entropy evolves in time and scales with 
the size of the system can be examined, and the physics of the growth of entanglement due to 
the propagation of quanta ought to be more interesting since there are more shapes available for 
region A. We examine the physics of the infinite strip and the round disc and indeed see rather 
interesting physics to compare and contrast, again following an evolution after a thermal quench 
modeled by a four dimensional Vaidya metric for gravitational collapse. A second issue we wished 
to address was how to do a different, non-thermal, type of quench in this AdS/CFT context, and 
to compare our results for it to the thermal case. We constructed a non-thermal quench by exciting 
electromagnetic sectors instead. There is a natural U{\) sector to which the physics couples and 
so starting with pure AdS4 we study and evolution of a shell (using a charged generalization of the 
Vaidya metric) which ultimately forms a zero temperature state at late times: an extremal black 
hole. So the quench is non-thermal, and the excitation in our new example is either by an external 
magnetic field or by a chemical potential for charge density, or an admixture of the two that can 
be freely dialed up using the four dimensional electric-magnetic duality of the four dimensional 
gravity system. 

In section 2 we give the general framework in which we are working, setting out our notation 
and the general form of the computation. We illustrate our methods by working out results for 
some familiar static geometries, re-deriving old results and adding some new ones (those for the 
case of extremal dyonic black holes). Section 3 presents our new results for the thermal quench in 
higher dimensions, where we notice several new features of the physics, including a kink in the time 
evolution of the entropy as it heads to saturation, for large enough system size. The case of the 
disc and the strip are thoroughly examined and contrasted. Section 4 presents the electromagnetic 
quench, discussing again many new features seen for the case of the disc and the strip. We conclude 
in section 5, and there is a brief appendix about the unitarity of the quench evolutions. 

2 General Framework 

In this section we present the setting for our computations of minimal surfaces in asymptotically 
AdS4 geometries. We will present some examples of static test geometries in section 2.3, with some 
familiar results and some new ones. This will allow us to set up notation and orient the discussion. 
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2.1 Gravity Background 

We consider the Einstein-Maxwell action: 

S = — ^ / d 4 xV^G (U-2A- F 2 ) + S ext , (3) 

1D7tG N J 

where the cosmological constant A = —3/R 2 sets a length scale R, and we will present the external 
sources we are interested in shortly. The equations of motion are given by: 

- ~G M „ (K-2A- F 2 ) - 2F llX F u x = 8irG 4 T^ , (4) 

-±=d» (y=GF^) = ^G^ c _ cxt) . (5) 
We consider a solution to these equations of the form [22J: 

ds 2 = -V(r,v)dv 2 + 2drdv + ^dx 2 , (6) 

R 

= r 2 ^ + qe{v? {v? ^ 
R* r r z 



Qe(v) i - _ </,„(> 

which correspond to having sources given by: 



4^G 4 J c M _ cxt = ^ V , torGtffi = (<7 e <fe + ? m g m - rro) . (9) 



This solution is a generalization of the Vaidya metric |231 [24"] . Note that in terms of the Hodge 
dual field, we can write: 

1 



d, V^G(*Fr = ^G 4 jy. , 4ttG 4 J* = ^ . (10) 



Also note that (m(v) , q e (v) , q m (v)) have dimensions of length. It is convenient to work in terms of 
a coordinate z: 

R 2 , \ 

* = — , n 
r 

such that we can write the metric as: 



ds 2 = ^ (-/(*, v)dv 2 - 2dzdv + dx 2 ) , f(z, v) = l- ^z 3 + qe(v) + J m{v) z 4 . (12) 



The time coordinate t of the dual field theory emerges (near the AdS boundary) as: 

v~t-z. (13) 
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2.2 Entanglement Entropy for the Strip and Disc 
2.2.1 The Strip 

We first consider calculating the entanglement entropy of an infinite strip, which we denote as region 
A in the dual field theory. We consider the area of a two-dimensional surface 7.4 that extends into 
the bulk that is joined to the boundary of A. We depict this in figure [ija). 

I 



< ► 




(a) The strip. (b) The disc. 

Figure 1: Diagrams of the two shapes we will consider for region A. 



We parameterize the surface with coordinates: 

e = x , e= y , (14) 

At z = e, very near the AdS boundary, (e is the UV cutoff of the dual field theory), the region A 
will be chosen such that: 

A-.= xe (--, -J , ye(0,L) . (15) 
with L taken to infinity. The area of this surface is given by: 

f i/2 r 

Area(7^) = L I dxyh , (16) 
J -11% 

where h a p is the induced metric on the surface given by: 

dx 11 dx u 

Kp-G, v — ^. (17) 

By the symmetry of the strip, we consider an ansatz of the following form for the embedding profile 
of j A : 

v = v(x), z = z(x), z(±£/2) = e, v(±£/2)=t-e , (18) 

such that the area is given by: 

. . , r „ 2 f £ / 2 , Jl - f(v, z)v'(x) 2 - 2v'(x)z'(x) , , 

Area( 7 ^ = LR 2 / dx^ y ' ) ' v 1 v 1 . (19) 

J -111 z ( x ) 
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We denote the integrand in equation (19) as the Lagrangian C of our system. Since there is no 
explicit x dependence, a constant of the motion is the Hamiltonian T~L (if we think of the coordinate x 
as a time coordinate), 

11 _ £ _ v >( x )_^l z '( x ) _ = 1 (20) 

'5v'(x) J 5z'(x) z (x) 2 s/l - f(v, z)v'(x) 2 - 2v'(x)z'(x) ' 

By the symmetry of the problem, and by choice of origin of x, we look for a turning point at x = 0, 
and at this turning point, we have: 

V '(Q) = 3/(0) = , 2(0) = , u(0) = v* . (21) 



Therefore, we have the result that H = z^ 2 . Using this result and equation (20), we can write the 
following conservation equation: 

1 - f(v, z)v\x) 2 - 2v\x)z\x) = . (22) 

By taking the x-derivative of this equation, and substituting for z"(x) in the two equations of 
motion, we get a single equation: 

- 4 + 4f(v, z)v'(x) 2 + 8v'(x)z'(x) + 2z{x)v" \x) - z(x)v'(x) 2 d z f(v, z) = . (23) 

Similarly, if we took the x-derivative of the conservation equation, and substituting for v"(x) in 
the two equations of motion, we get a single equation: 

4f(v, z) 2 v'(x) 2 - f(v, z) (4 - 8v'(x)z'(x) + z(x)v'(x) 2 d z f(v, z)) 

-z(x) (2z"{x) + v'(x) (2z'(x)d z f(v, z) + v'(x)d v f(v, z))) = . { ' 

The "on-shell" area is now given simply by: 

r£/2 z 2 

Area(7^) = 2LR 2 / dx— ^ . (25) 
Jo z \ x ) 

2.2.2 The Disc 

We may now consider the disc. We consider a disc of radius I, and take the two dimensional surface 
7^4, which we depict in figure [jjb), to be parameterized by: 

e=r, e = 4>- (26) 

By the rotational symmetry of the problem, we can consider an ansatz of the form: 

z = z(r) , v = v(r) . (27) 
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This gives for the area: 

Area(7^) = 2ttR 2 / dr-^yjl - f(v, z)v'(r) 2 - 2v'{r)z'{r) . (28) 
Jo z \ r ) 

Because there is an explicit r appearing in the Lagrangian, we do not have a conserved quantity 
associated with the disc as we did for the strip. However, by the symmetry of the problem we again 
must have: 

Z '(Q) = „'(o) = . (29) 



The general equations of motion derived from the extremizing of the area in equation (|28j) are 
straightforward to derive. We do not write them here, in general, as they are rather long. Later, 
specific cases will be written. 

2.3 Static Examples 

Following are computations for the entanglement entropy for three static spacetime examples. In 
each case we present the results for the strip and the disc. 

2.3.1 Pure AdS 4 

Let us calculate the entanglement entropy for pure AdS 4 [7]. This means taking f(v,z) = 1. 
Although it is simpler to exchange v for t, let us proceed with the given choice of coordinates. A 



solution to the equations of motion (23) and (24) is simply to take: 



v '(x) = - z '{x) . (30) 



Then the conservation equation (22) allows us to write: 

dz 



* ± V?- 1 - (31) 

where the plus sign is taken for x < and the minus sign is taken for x > 0. We can use this result 



to find a relationship between z* and t by simply integrating equation (31): 

4 \ -1/2 

2 Je - 

Changing coordinates to u = z A /zf, we can write this integral in a more convenient form: 



fU dz(^-l) • (32) 



2 ~ T 



(33) 
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where B(a,b) is the Euler beta function, and B(x;a,b) is the incomplete beta function. In this 
expression, we can take e — > without encountering any divergences, so our result is simply: 

£ _ VHW4) 

2~ z * r(i/4) • l<34j 



Using the same procedures, we can write our "on-shell" area as: 

4 , , 2LR 2 f z * , z, 2 / z 4 \ 

Area( 7 ^) = / ^ 4 1 - - . (35) 



Changing coordinates to u = z /z%, we can write this integral in a convenient form: 

/.'a.-- ^<i-.r*-^(s(-i.i)-*(£4i)) • (36) 

Expanding our expression for small e (in order to capture the divergent contribution from our UV 
cutoff), we get our final result for the area [7]: 

Areaf^l - ^ + * L -2R L ( R + R ^(-l/4)r(3/4) \ 

Area( 7 ^) - -— + 2 ^ r(1/4) - 2«*< + ^ r(1/4)2 j • (37) 

The divergent first term is proportional to the boundary of A. We will be generally interested in 
the finite contribution to the entropy given by the second term, which we denote s: 

1 /. . . 2R 2 L\ . , 

Area( 7 ^) (38) 



4G N 

We will work in terms of the rescaled finite entropy s, a function of terms of the dimensionless 
length £ = £/R, 

~ 4Gn (M 

S= 2RL S > (39) 

which we plot in figure [2] 

We can proceed in a similar fashion for the disc. A solution to the equations of motion is 



to take the result of equation (30), which then gives us one equation for z(r): 

2r (1 + z'{r) 2 ) + z(r) (z'(r) + z'(r) 3 + rz"{r)) = . (40) 
A solution to this equation is: 

z(r) 2 + r 2 = £ 2 , (41) 

which describes the extremal surface as a half-sphere (of radius £) in the bulk. The area is given 
by: 

r-t 

Area(7^) = 2irR 2 £ / dr— T — . (42) 
Jo z[rY 
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Figure 2: The finite contribution to the entanglement entropy as a function of I = l/R in the pure AdS4 
background. 

We use the fact that dz = —rdr/z and z* = £, and we change coordinates to y = zji to write the 
integral as: 

Area( 7 ^) = 2irR 2 [ % = 2itR 2 ( - - 1 ] . (43) 

If we define the finite contribution to the entropy by: 

If. . , 27ri? 2 A 

s= 4^r ea( ^ — ■ (44) 

and consider the rescaled entropy s 

then we find that s = — 1. Note in particular that there is no dependence on £, the size of the 
disc [7]- We emphasize that this is a significant difference from the strip and is a simple result of 
scale invariance in the theory when a shape that preserves conformal invariance on the boundary 
is taken for A. 

2.3.2 AdS 4 Schwarzschild 

To study the case of AdS4-Schwarzschild[7], we take f(v, z) = 1 — 2mz 3 /R 4 . In this geometry, the 
event horizon is located at a radial distance zq given by: 

We once again start with the strip. A solution to the equations of motion is to take: 

v(x) = t + g(z(x)) , d z(x) g(z(x)) = -/(z(x))- 1 . (47) 
The remaining equation of motion is given by: 

- 2z(x)z"(x) + z'{xf (-4 + zix) 9 ^^^ ) - Af(z) = . (48) 
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The conservation equation is now given by: 



z'(x) 2 

JW))~W) 4 



1 + 



which allows us to write: 



dz 
dx 



±a//(*(*)) 



z(x) 4 



In turn, the area can be written as: 

Area(7^) 
and the length is given by: 



2LR 



2 rz 



dz^ f(z) 1 



-1/2 



(49) 



(50) 



(51) 



dz-, f(z) 1 



-1/2 



It is convenient at this point to define dimensionless variables: 



z x 
z = — , X = — 
ZQ ZQ 



such that: 



and 



Area(7^) 



2LR 



2 rz 



dzz 2 



(l-^)U 



-1/2 



z, ~2 



dz- 



(l-z 3 )(l 



~ 4x N -1/2 
Z ■ 



(52) 



(53) 



(54) 



(55) 



Since there is no analytic expression for the area, we wish to express the finite contribution to the 



entanglement in a way that is consistent with our pure AdS4 result using equation (38). Further- 



more, we can define a rescaled entropy s as in equation (|39|): 

. R 



dzz 2 



W) i 



-1/2 



(56) 



We present a plot of this finite contribution to the area in figure [3j In particular, note that for large 
£, the result for zqs/ R scales as £/2 which is exactly what is needed to recover the entropy-area 
law, i.e.: 

(57) 



R 2 

(4Gn)s ~ ~2^£ , for £/zq large . 



We can perform a similar computation for the disc. The area can be written as: 



Area(7A) = 2ttR 2 



r I z'(f) 2 
df—\ll + 



1-z 3 " 



(58) 
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zo. 
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R 




Figure 3: The finite contribution to the entanglement entropy for the strip as a function of £ = 1/zq in the 
AdS4-Schwarzschild background. 



We write the finite contribution to the entropy s as in equation (44), and we write the rescaled 



entropy s as in equation ( 45 ) : 



r , z'(f) 2 
dr—\/l + 



(59) 



/o z* V 1 — z 3 e 

We present the behavior of the area in figure]^} We see that for small I, the contribution matches 
that of pure AdS4, as the surface remains close to the AdS boundary. For large £, the behavior of 
the area is quadratic in £ as is to be expected, since most of the contribution will come from the 
surface behaving as a disc at the event horizon. Again we find that for large £, the entropy s scales 




Figure 4: The finite contribution to the entanglement entropy for the disc as a function of £/zq in the 
AdS4-Schwarzschild background. 

as £ 2 /2, which again recovers the entropy-area law: 



(4G N )s 



R 2 



ir£ 2 , for l/zo large . 



(60) 
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2.3.3 Extremal AdS 4 Dyonic 



We now consider the extremal AdS4 dyonic black hole. The entanglement entropy for this static 
case has not been studied (as far as we know) in the literature, and so the results are new. This is 
a useful preparatory computation for the electromagnetic quench we will perform in a later section. 
We have: 

f(z) = l-4 



where we have used that: 



Q 



2,2 
Qe ' Qm 



—R 2 
2z 



2 



2R 4 



in 



(61) 



(62) 



Equations (47), (48), (49), (50), (51), and (52) remain unchanged. We can define the finite part 



of the entropy using equation (38) and the rescaled entropy using equation (39). We present the 



finite contribution to the area in figure [5} As in the case of the Schwarzschild black hole, for large £ 
we recover the entropy-area law with a linear dependence in I. Performing a similar calculation 




Figure 5: The finite contribution to the entanglement entropy for the disc as a function of £/R in the A0IS4 
Dyonic black hole background. 



for the disc gives the results in figure |6| We see again that for I large, the behavior is such that we 
recover the entropy-area law, with a quadratic dependence in i. 

A new feature of this background is that the area integral exhibits a logarithmic divergence 
as z+ approaches the event horizon. This is a simple result of the fact that the event horizon has a 
double zero. This feature will play a role in the evolution of our electromagnetic quench discussed 
later. To see this more explicitly, as an extreme example, consider the embedding with z* = zq. 
This embedding consists of a piece that extends along the event horizon, and then two pieces that 
extend from the AdS boundary at x = ±£/2 (for the strip) to the event horizon. We present a 
depiction of this embedding in figure [7| The contribution to the area of the pieces that fall towards 
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Figure 6: The finite contribution to the entanglement entropy for the disc as a function of t/R in the A0IS4 
Dyonic black hole background. 

£ 




Z — Z* 




z = e 



z = z 

(a) Embedding with z, < zq (b) Embedding with z* = zo 

Figure 7: Diagrams of the two types of embeddings in the presence of a black hole. 



Z = Zq 



the event horizon is given by: 



Area( 7 ,0 = 2L / -^f(z) 



-1/2 



Near the event horizon, we can write: 



f(z -> Zq) = (ZQ ~ zf , 



and so we have: 



Area(7^) 



Z( » dz 



ZQ- Z 

which contributes a logarithmic divergence to the area. 



(63) 

(64) 
(65) 



3 The Thermal Quench 

We begin by extending the very interesting work of refs. 0120] to one dimension higher, performing 
a thermal quench. We will get several new results, and contrast them with several more new results 
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that we obtain for a new, electromagnetic, holographic quench in a later section. The idea of the 
holographic quench studied in refs.|8, 20] (see also ref.[25]) is that the system is simply AdS at 
early times, and then at some later time it is strongly perturbed by a shell of collapsing null dust, 
which subsequently forms a Schwarzschild black hole at late times, representing a thermal state 
in the holographically dual theory [2T]. The point of the exercise is to study the evolution of the 
entanglement entropy over time after the quench^ Generically, it saturates after a characteristic 
time. In 1+1 dimensions the behaviour is rather simple, with a linear rise in the entropy followed 
by a distinct leveling off to saturation, as observed in a magnetically quenched system in ref. |3j, 
and a thermally quenched system using holographic duality to asymptotically AdS3 geometries in 
refs.[U [20]. The slope of the line reflects the fact that the effects of the quench propagate at the 
speed of light in the single dimension available. We shall see more complicated possibilities in 
higher dimensions and in the next section when we do an electromagnetic (non-thermal) quench. 
To achieve the quench, we consider: 



f(v,z) 



1 



zo(v) 



m[v) 



zo(v)~ 



(66) 



and take the following functional form for zq(v) 

. 1 + tanh(v/uo)\ 

zq(v) - - 



(67) 



This form for zq{v) is such that the background is pure AdS4 in the infinite past and an AdS4~ 
Schwarzschild black hole in the infinite future. See figure [8} The parameter z^ can be thought of 



Zo(v)~ 



1.0 




0.8 




0.6 




0.4 




fi.2 





Figure 8: The shape of the mass function m(v), generating the quench from pure AdS4 to a finite mass 
black hole. 



1 This evolution is indeed unitary, as it ought to be [20] ■ See the appendix for some clarifying comments on that 
issue. We thank E. Lopez for a discussion on this point. 
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as being related to the final temperature of the background at late times: 

Air 00 ' 



(68) 



and the parameter vq is related to the time over which the quench occurs. As in the case of the static 
AdS4-Schwarzschild black hole, we rescale z and x by z^ to define dimensionless coordinates z 



and x. We proceed to solve equations (23) and (24) numerically using a shooting method with the 
following initial conditions: 



ff(Q) = z'(0) = , 5(0) = 5* , v(0) = v* . 



(69) 



Near the AdS boundary (i.e. when x = £/2 for the strip or when f = £ for the disc), we extract 
the following information: 

~z(£/2) = ~e, v(£/2)=t-e, (70) 
Because of the form we choose for our quench, the quench starts at approximately t ~ — 2vq. 

3.1 The Strip 



After the rescaling, the area of the surface given in equation (25) is given by: 
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Area(7^) 



2LR Z 



dx 







~z(x) A 



(71) 



We define the finite contribution to the entanglement entropy s as in equation (38), but we modify 



our definition of the rescaled entropy from equation ( 39 ) by subtracting the pure AdS4 result we 
have found in section [2,3.11 



R 

Zoo 



dx 



z(x) A 



tt r(-i/4)r(3/4) 



2£ 



T(l/4f 



(72) 



and with this subtraction we have chosen that the entanglement entropy s starts at zero in the 
infinite past. We present some results in figure [9| showing the evolution for a range of quench 
speeds, set by %. We see that the curves almost perfectly overlap when vq is small enough, except 
in the region before saturation is reached, where the curve with smaller vq appears to be sharper. 
The invariance of the early evolution suggests that regardless of how fast the quench is done, the 
early evolution is dictated by a time constant that is only dependent on the details of the plasma 
described by the dual 2+1 dimensional field theory. This behavior was also observed in refs.[8l 120] 
for one dimension fewer. Notice that for large vq, the evolution of the entropy appears to follow 
the evolution of the quench rather closely. For short quenches the system relaxes to its new state 
more slowly. 
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Figure 9: The evolution of s for different quench speeds for I = 1. 



Exploring the behavior as we change the strip width £, we note several things. First, 
the early evolution appears to be independent of (or only weakly dependent on) I. This can be 
understood from the geometry of the strip. As £ changes, the area at the boundary of A does 
not change, so the initial propagation of excitations from region A to region B, contributing to 
the entanglement, is not affected by the strip width. Second, we see the appearance of a new 



phenomenon, the appearance of a swallow tail. See figure 10 The swallow tail is a feature of 
our numerical method for searching for minimal area solutions with a given UV and IR boundary 
condition at a given time t. The method reveals the presence of multiple (at times up to three) 



solutions to the problem for a given time. By studying figure 11(a), we see that the multiplicity 
corresponds to the presence of solutions that have a below and above the apparent horizon. The 
solution is always the one with the lowest area, and therefore the transition leads to a kink in the 



entropy dependence on time. Figure 11(b) explicitly shows the three embeddings at the transition; 
the transition is from an embedding with several sharp features to a smooth one. 

This kink in the evolution of the entanglement entropy is a striking new feature in time 
evolution. Notice that the branch of solutions joined on to at late times are ones for which the 
turning point of the minimal surfaces, £*, is well above the final horizon, Zqo, giving an s which is 
essentially the saturated late time value. In other words, the entropy reaches saturation and then 
abruptly flattens out for large enough strip size £. For smaller strip size, the transition is also fast, 
but not so fast as to have a discontinuity in the derivative. Since the rate of approach to saturation 
is not strongly dependent on the region size £, while the saturation value itself depends on £, it is 
reasonable to expect that upon nearing saturation, the evolution turns over and flattens out toward 
the saturation value. In cases where the saturation value increases fast enough with £ (linearly with 
I for the strip) while the rate of approach remains constant or weakly dependent on £, this turnover 
will be more pronounced for large enough £. In these large N systems we are studying, this is, we 
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presume, the role of the kink, which as we have seen appears for large enough t. We expect that 
the sharp transition would be softened by 1 /N corrections upon leaving the gravity approximation. 

Our results suggest that the picture proposed in ref. [3], which arose in the 1+1 dimensional 
case, needs modification here. In their picture, the entropy is discussed in terms of left and right 
moving pairs in the field theory. Saturation occurs once all pairs have had enough time to have 
traveled a distance of order the size of the interval. For our higher dimensional geometries we see 
saturation before all the analogous pairs are separated enough (here for the strip, for example, 
there are directions for which pairs can travel for extremely long times before leaving the strip). 
This suggests that the entanglement entropy is not built from simple pairwise contributions for 
the higher dimensional cases under consideration here. The bulk of the contributions from those 
longer-traveling quanta must have been already accounted for, and the system nears saturation 
without them having departed the starting region. 
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Figure 10: Results for thermal quench for vq = 0.01. 



Turning to the intermediate time region where the entanglement entropy grows rapidly, we 
show in figure 12 ) the behavior of the saturation time for various I. Near i = 0, the best fit line has 
slope 0.77, yet for much larger £, the slope is closer to 0.65 (Note that there is some flexibility in 
this number depending on how one defines the edges of the growth period to extract the saturation 
time t sa t)- Interestingly, therefore, it seems that the saturation rate for the thermal case grows 



faster for larger £. This is traceable to the early time behaviour of the curves. See figure 10 For 
small I («.e.< 2), saturation is reached well before we enter the linear phase , so we may expect 
differences in the saturation time for those cases versus large I. 

It is important to note that our saturation time is a departure from the result of 1/2 for 
the 1+1 CFT result of ref. [3], reproduced in the holographic thermal quench studies of refs.[8l [20]. 
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(a) Evolution of 5* 



0.5 1.0 1.5 2.0 % 

(b) Embeddings at the transition 



Figure 11: Behavior of z* for £ = 4 and vo — 0.01. In (a), the red dashed line is where the transition occurs 
and the apparent horizon is at z = 1. In (b), we display the two embeddings at the transition, where the red, 
dot-dashed curve corresponds to the embedding with z* = 2.0669, the green, dashed curve corresponds to 
the embedding with z* = 1.2010, and the blue, solid curve corresponds to the embedding with z* = 0.9985 

Note that here the excitations do not necessarily travel a distance of order 21 during this time. 
Because of the geometry of the strip, some excitations can travel much greater lengths. Therefore, 
the overall effective speed of saturation of s should be expected to be less than the speed of light, 
as we have seen in our results for the strip. 



? sat -2v 




Figure 12: The saturation time as a function of £ for thermal quench of the strip. The solve of the best fit 
line for small £ depicted with the solid blue line is 0.77. The slope of the best fit line for large £ depicted 
with the red dashed line is 0.65. 
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3.2 The Disc 



We can perform similar calculations for the disc. The area of the two dimensional surface is given 



by equation (28), and we define: 

df— 



■y/l - f(v,z)v'(r) 2 - 2v'(r)z'(r) 



+ 1 



(73) 

We present some results in figure [13} There are several differences from the strip case. First, we 
see that the early evolution strongly depends on £. Once again, this can be understood from the 
geometry of A. The region near the boundary of A has a size that grows with I and so the early 
evolution of the entanglement should reflect this dependence on geometry. Another key difference 
from the strip case is that no matter how large t is, we do not see the kink in the transition 
to saturation. In line with our reasoning from the previous section, the saturation value of the 
entanglement entropy and the rate of approach to saturation increase in lock with I (unlike for 
the strip where the rate was independent of €) allowing for a smooth evolution towards saturation. 
Furthermore, we see that the saturation time grows linearly (with a slope of unity) with I. This fits 
with the fact that the disc's finite extent in all directions means that there is a maximum travel time 
for any excitations as a result of the quench to leave A and contribute to the entanglement entropy. 
So in this way we recover that the rate of saturation of the entanglement entropy is the speed of 
light, as was observed in lower dimensions. In this sense, the disc is a more natural generalization 
of the 1+1 dimensional case than the strip. 
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(a) Thermal quench, with vq = 0.01. 
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(b) Saturation time for varying t, with vq — 0.01. 



Figure 13: Results for thermal quench with a disk. The slope of the saturation time is almost exactly unity. 



20 



4 The Electromagnetic Quench 



It is possible to cleanly study an entirely different type of quench in this system. This non-thermal 
quench is achieved by turning on electric and magnetic sources of the theory. We choose the system 
to be entirely non-thermal by approaching the extremal (i.e., zero temperature) black hole solution 
at late time^j The electric components of the background at late times can be read off as non-zero 
charge density and chemical potential in the dual 2+1 dimensional field theory [261 127]. whilst the 
magnetic components represent a background magnetic field |28j. In fact, electric-magnetic duality 
invariance of our system means that different choices of electric or magnetic sources can be rotated 
into each other, and we will often think of the system as undergoing a purely magnetic quench 
for convenience, but admixtures of electric (chemical potential) and magnetic (background field) 
quenches are captured by all of our results. 

The collapsing charged geometry has a metric similar to the extremal AdS4 metric from 
the section [2. 3 .3[ but with the function / now taken to be: 



3 / \ 4 



f(v,z) = l-4 — +3 — - , (74) 



taking: 

q e (v) 2 + q m (v) 2 = 4 , m(v) = . (75) 



3R 6 , . 2R 4 

z {vp zq[vY 



We use the same quench profile for zq(v) as that given in equation (67). Note that here z^ sets 
the final strength of the magnetic field (in the duality frame where we choose everything to be 
magnetic) . 

4.1 The Strip 

We proceed in a similar fashion as in the case of the thermal quench. The area is given by 



equation (71), and we define the finite entanglement entropy s by equation (72). Let us consider 
first the effect of quench speed, set by vq. We present an example in figure 14 As with the 
thermal quench, it seems the general behavior of the entanglement entropy for small enough t is 
qualitatively universal. When v$ is small enough, it appears that the rise time to linear behavior 
does not depend on vq, suggesting that there is an intrinsic time scale for the dual plasma to react 
to changes in the medium, as in the case of one dimension fewer [2D]. It is interesting to note, 
however, that the time scale in the current system appears to be much smaller than that of ref.[20j. 
As in the thermal case, for large vq, the evolution of the entropy appears to follow the evolution of 



2 Since the extremal solution is a highly degenerate point in the space of black hole solutions, it is clear that 
forming one by collapse is a delicate process. However, it is possible, at least formally. It turns out that the novel 
physics that we observe in this system does not appear to depend upon the details of the process of formation. 
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Figure 14: The entropy and magnetic field for £ = 1. 



the quench rather closely. We consider the effect of strip size £, presenting some results in figure 15 
The early time evolution is independent of £ just as we saw for the thermal case. Furthermore, for 
large enough £, there is again the appearance of a swallow tail, controlling the rapid flattening out 
to saturation. The resulting kink in the entanglement entropy has a similar origin to that discussed 



in section 3.1 although in this case the situation is more severe since the rate of approach is slowed 
down by the logarithmic behavior. An increase in £ give an increases in size for the swallow tail, 
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Figure 15: The entropy for various values of £ and vq = 0.01. 
and eventually the transition occurs even before the beginning of the swallow tail region is ever 



reached. See an example in figure 16 



A new important feature is a departure from the linear behavior of the evolution to a 
logarithmic growth. This behavior is inherited from the logarithmic behaviour we saw in the static 



extremal case in section 2.3.3, which has its origin in the double zero at the horizon. On the dual 
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field theory side, this means significantly new behaviour for the excitations as they propagate after 
the quench. Presumably, the charged excitations are slowed down by electromagnetic interactions, 
either from screening in the presence of non-zero charge density or due to propagation in the 
background magnetic fielcJ^J 




Figure 16: The entropy for I = 5 and vq = 0.1. The intermediate time behavior of the area is well 
approximated by Int. The dashed line corresponds to the second embedding to which the system transitions, 
showing entanglement saturation. 



We can study how the saturation time i sat depends on I, We present these results in 
We find that the slope of the line near t = is approximately 0.77 (note again that the 
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figure 

slope depends slightly on how one decides to calculate t sa t because of the roundedness of the solution 
before it saturates), which matches our earlier result for the thermal quench. Upon reaching the 
point where the jump between embeddings starts to occur, the growth in t sat appears to be well 
fitted with an exponential from that point onwards (to be expected from the logarithmic growth 
of the area with time that we commented on earlier). 

4.2 The Disc 

We proceed in a similar fashion as in the case of the thermal quench. The area is given by 



equation ( 28 ) , and we define the finite entanglement entropy s by equation ( 73 ) . We present some 



results in figure 
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Again, the early evolution has strong I dependence in a manner attributable to 
the disc geometry. Also present is the logarithmic evolution, arising from the near horizon region 
of the geometry. Furthermore, for large enough £, there is an appearance of the swallow tail. Recall 



that this was absent for the thermal quench for the disc. In line with the discussion in section 3.1 



3 Logarithm growth of entanglement entropy has been observed in (l+l)-dimensional systems quenched with an 
external field (see e.g. ref . 129j ) . It would be interesting to explore the possibility of a connection to our results. 
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Figure 17: The saturation time as a function of I for vq = 0.01. The slope near £ = is approximately 0.77. 

its appearance here is due to the slowing (as opposed to increasing in the thermal case) rate of 
approach to saturation with increasing t due to the logarithmic behavior. 

5 Conclusions 

We have studied two types of quench for the (2+l)-dimensional theory dual to AdS4, computing 
the evolution of the entanglement entropy for an infinite strip region and a round disc region. Our 
quenches were a thermal quench dual to a shell of uncharged matter (which forms a Schwarzschild 
black hole) and a non-thermal quench, formed from a shell of charged matter (which we tuned to 
form an extremal (i.e., zero temperature) dyonic black hole. 

We observed several interesting phenomena. The results for the disc, for the thermal quench, 
provided the most direct analogue of earlier observed results from (1+1) dimensions [5} IBj 120] : There 
was linear growth in time of the entanglement entropy on its way to saturation, with a slope 
consistent with the interpretation that the quanta responsible for entanglement were propagating 
with the speed of light. The strip deviated from this in two ways. The first, that the effective 
growth rate (while still linear) was slower than the disc's was naturally attributable to the semi- 
infinite geometry of the strip. The second was more surprising. For large enough strip size there 
was a kink in the entropy-time curve corresponding to a rapid turnover at saturation. The source 
of this on the geometry side is quite natural (it is due to the availability of multiple branches of 
candidate minimal surface solutions, forming a swallow tail shape — there is a jump from one 
branch to another at late enough times), but it would be interesting to understand its origins on 
the field theory side. The kink suggests a modification of the picture proposed in ref. [3], whereby 
not just entangled pairs contribute to the entanglement entropy. (Note that it is probable that 
such a kink gets softened beyond the supergravity approximation.) 
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(a) Magnetic quench vq = 0.01. 



(b) Swallow tail, £ = 2.75, v = 0.01. 
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Figure 18: Results for magnetic quench with a disc. The slope of the saturation time for small £ is almost 
exactly unity. 



The kink is also present for both the disc and the strip for the electromagnetic quench. One 
of the features of our novel quench is that the entanglement entropy grows only logarithmically 
with time, a feature that may be due (on the field theory side) to screening effects for the charged 
quanta. It is interesting that kinks appear when the effective propagation speed contributing to 
the entanglement entropy growth is subluminal for both thermal and electromagnetic cases, being 
absent for the thermal quench disc case. We do not know if this is a coincidence, but are encouraged 
to pursue further study of this issue. 

In interpreting our results, we have assumed that the prescription to calculate the entan- 
glement entropy is correct. New features like the kink have been discussed from this point of view 



(see section 3.1). We recognize that, until a better understanding of the prescription is obtained, 
it is possible that the kink may be interpreted as pathological, indicating a breakdown of the 
prescription to calculate the entanglement entropy for rapid quenches and large system size. 
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The logarithmic growth is again something that would be interesting to understand further 
on the field theory side. On the gravity side it comes from surfaces that probe near the extremal 
horizon, which has a double zero in the metric function. A throat opens up at the horizon, producing 
a smooth geometry that is AdS2 x M 2 . It is tempting to speculate that for our electromagnetic 
quench, the growth of entropy in the intermediate region might be captured by an effective lower 
dimensional theory captured by AdS2 physics, as has been shown to happen for aspects of the 
physics of studying Fermi surfaces holographically in a related context [301 13"T]. 
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Appendix: A Note on Unitarity of Quench Evolution 

In order to show that the evolution of the quench and subsequent relaxation process is unitary, we 
would like to show that throughout the evolution the following condition is satisfied: 

S A = S B . (76) 

It is useful to start with some static examples before considering an evolving background. Let us 
begin with the case of pure AdS4. To compute there are two possible surfaces to consider [12J, 



which we depict in figure 19 The embedding in figure 19(a) can be shown to be divergent by first 



considering it to be a sum of two finite intervals and then sending the outer edges of the intervals 



to infinity. The embedding in figure 19(b) has a piece that gives S_a and the contribution from the 
piece at z = oo is zero. Therefore, we find that AdS4 corresponds to a pure state. 

Let us now consider a static Schwarzschild black hole. Repeating the calculation described 



above, we find that the second piece of figure 19(b) must lie not at z = oo but at z = Zq, the 



horizon. Because of this fact, its contribution gives a non-zero value; in fact it gives an entropy of: 

1 ( R2 f j j A a Area(E.H.) 



SB = SA+ 4G^ [4 J dxdV ) = SA+ 4G n ' ' (7?) 

where Area(E.H.) is the area of the event horizon. Therefore, we find that Sb 7^ and so the 
system is not in a pure state. Therefore, any static solution with an event horizon corresponds to 
a mixed state (see e.g. ref. [32 1). 
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Figure 19: Diagrams of the two types of embeddings to calculate entanglement entropy of B. 



So let us now consider our evolving metric. If we consider a solution that corresponds to 
the embedding of figure |19(b)| a solution to the equations of motion for the second disjoint piece 
is to simply have: 

z(x) = oo , v(x) = — oo . (78) 

We note that for v(x) = — oo, our metric does not have a singularity at z{x) = oo, so there is no 
danger for the solution to lie there. In particular, this solution is (in a sense) simply the solution 
that we would get for AdS4 and gives zero contribution to the entropy. The first piece will give us 
the entropy of A. Therefore, for the evolving case, we find that we have: 



(79) 



Note that this only works because we are assuming the second piece remains static in the sense 
that it remains at v = — oo whereas the other piece evolves. Therefore, with this in mind, we find 
that the evolution is unitary. 
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